NASA TR R-41 


NASA TR R-41 


NATIONAL AERONAUTICS AND 
SPACE ADMINISTRATION 

TECHNICAL REPORT 
R-41 

NORMAL COMPONENT OF INDUCED VELOCITY FOR 
ENTIRE FIELD OF A UNIFORMLY LOADED LIFTING 
ROTOR WITH HIGHLY SWEPT WAKE AS 
DETERMINED BY ELECTROMAGNETIC ANALOG 

By WALTER CASTLES, JR., HOWARD L. DURHAM, JR., 

*nd JIRAIR KEVORKIAN 


1959 


For sale by the Superintendent of Documents, U.S. Government Printing Office, Washington 25, D.C. Yearly subscription, $10; foreign, $13; 

single copy price varies according to size --------- Price 40 cents 




TECHNICAL REPORT R-41 


NORMAL COMPONENT OF INDUCED VELOCITY FOR 
ENTIRE FIELD OF A UNIFORMLY LOADED LIFTING 
ROTOR WITH HIGHLY SWEPT WAKE AS 
DETERMINED BY ELECTROMAGNETIC ANALOG 

By WALTER CASTLES, JR., HOWARD L. DURHAM, JR., 
and JIRAIR KEVORKIAN 


Georgia Institute of Technology 


i 




TECHNICAL REPORT R-41 


NORMAL COMPONENT OF INDUCED VELOCITY FOR ENTIRE FIELD OF A 
UNIFORMLY LOADED LIFTING ROTOR WITH HIGHLY SWEPT WAKE AS 
DETERMINED BY ELECTROMAGNETIC ANALOG 1 

By Walter Castles, Jr., Howard L. Durham, Jr., and Jirair Kevorkian 


SUMMARY 

Value# of the nondime ns ional normal component 
of induced velocity throughout the flow field of a 
uniformly loaded lifting rotor operating in the upper 
half of the helicopter speed range are presented in 
the form of graphs and tables. The tabulated data 
are for rectangular grids of points located in azimuth 
planes situated at 80° increments of azimuth angle. 
The grids extend a distance of 4 rotor radii in both 
the vertical and radial directions. Values at points 
in the rotor plane were computed by means of the 
Biot-Savart relation using the assumption that the. 
wake-vortex distribution consisted of a uniform , sem i- 
infinite, elliptic cylinder. Values at points not in 
the rotor plane wen obtained experimentally by 
'measurements oj the field strength about an electro- 
magnetic-analog model of the wake-vortex system. 

( Jomparisons of computed and experimental analog 
values for the normal component oj induced, velocity 
both in the plane of the rotor and in the lateral plane 
perpendicular to the rotor plane are presented. The 
agreement between the computed and experimental 
analog values indicates that the latter are sufficiently 
accurate for engineering purposes. 

The results should be useful for estimating the 
induced velocity distribution about lifting rotors in 
general and for synthesizing the distributions over 
the rotor disk for the case of any specified nonun if or m 
hading. 

INTRODUCTION 

In order to determine the performance and air 
load distribution of a lifting rotor, it is necessary 
to know the induced flow field in the vicinity of 
the rotor, the component of velocity normal to 
the plane of the rotor being of particular interest. 

1 SuiHirsodcs NT A C A Technical Note I2MK by Walter Castles, Jr., Howard I 


To make rotor-flow-field computations mathe- 
matically tractable, it is usual to approximate the 
actual wake-vortex system by one having regular 
geometric properties. In general, however, for 
even the simplest of wake geometries the calcula- 
tions are tedious and prohibitively lengthy unless 
high-speed computing facilities are available. 
Alternatively, there is an approach to the problem 
making use of the perfect analogy between the 
induced How field associated with a vortex filament 
in a perfect fluid and the magnetic, field in space 
associated with a current-carrying wire. Thus it 
is possible to construct an electromagnetic analog 
in the form of a wire model of a given vortex (•on- 
figuration. Point measurements of magnetic- 
field strength in the associated magnetic field then 
afford a description of the analogous induced 
velocity in the fluid velocity field, as shown in 
reference 1 . 

The principal objective of the present paper is 
to present in the form of tables and graphs the 
experimental values for the tiondimensional normal 
component of induced velocity which were ob- 
tained by means of an electromagnetic-analog 
model of the wake from a rotor operating in the 
upper half of the flight speed range. The method 
employed was in many respects similar to the 
procedures described in references l and 2. Sur- 
veys were made of the normal component of 
induced velocity in several azimuth planes per- 
pendicular to the plane of the rotor beginning with 
the longitudinal plane of symmetry and proceeding 
in B()° increments of azimuth angle. 

Another objective is to supplement and extend 
the results of references B and 4 by presenting 

Durham, Jr., and Jirair Kevorkian 1!U>K 
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additional computed values, obtained by means 
of a digital computer, of the normal component 
of induced velocity in tin* rotor plane. This pro- 
gram was carried out along with the magnetic- 
analog measurements and atlorded reliable cheek 
points for comparison ot results. The computed 
data furnished values for the induced velocity at 
space points located such that physical inter- 
ference between the piekup coil and wake model 
prevented held measurements and also at points 
near the model coils where the gradient of the 
local magnetic field was large. 

The analysis presented herein concerns the flow 
field associated with a uniformly loaded lifting 
lot or and assumes that the wake-vortex system 
has the form of a uniform, semi-infinite, elliptic 
cylinder composed of a very large number of 
circular vortex ring elements arranged in such a 
way that the circulation per unit length of the 
vortex sheet is constant. This assumption implies 
that the induced How associated with the vortex 
system is a potential flow and, as such, has a 
perfect- magnetic analogy as pointed out by 
equations (2) and (d) of reference 1. 

This investigation was conducted at the Georgia, 
Institute of Technology under the sponsorship 
and with the financial assistance of the National 
Advisory Committee for Aeronautics. 

SYMBOLS 



coefficient of cosine term in Fourier 
scries for blade flapping angle 

MR 

output meter reading, dh 

in 

tangent of wake angle x 

V 

any ]>e>int Y{X Oi Y u ,Z 0 ) in rotor flow 
iielel 

R 

rotor radius 

Ro 

radius e>l' point P from Z— or rotor 
axis 

Rc,i/-c 

>'o RJlt 

cylindrical coordinates of a point on 
the curve of intersection of plane 
^-Constant with wake-vortex 
cylinder 

r 

velocity of helicopter along llight 
path 

r, 

normal component of induct'd ve- 
locity at P 

V 

normal component of induced ve- 
locity at center of rotor plant* 


X,Y,Z 

A \,Y ut 2 0 
z 0 ZJii 

a 

at) 


coordinates of a wake-vortex sheet 
element as measured relative to the 
t ip-pat h-phme axes 
coordinates of point. P in rotor flow 
field 


angle of attack of plane of zero 
feathering 

angle of attack of rotor plane, 


a — Hi 

6 azimuth angle of wake-vortex sheet 

element having length <ls meas- 
ured from negative X- or up- 
wind direction 
\ P = (V sin a 9 —o)/&R 
/* P =(V :*os a v )/QR 

X wake angle measured between nega- 

tive Z— or rotor axis and wake 


axis 

i p azimuth angle of point P measured 

from positive X— or downwind 
direction 

£2 angular velocity of rotor blade's 

Subscripts: 

(! curve of intersection formed by plane* 

^ = Constant with wake-vortex 
cylinder 

N search-coil normalizing point for 

which computed velocity ratio was 
known 

Y point P 


THEORETICAL ANALYSIS 


Under the assumption that the wake- vortex 
distribution takes the form of a uniform, semi- 
infinite, elliptic cylinder, it was shown in reference* 
4 that the ratio of the* normal component of in- 
duced elocity at any point- P to Unit at the center 
of the color is given by 



2tt Jo n 7j (yV/— />) 


(16 


0 ) 


where lie wake geometry is given in figure 1 and 
A -1 +r w cos ($—6) 

B m cos fl/y 1 -j -v/r 
C T =l+r tf a -l"3* 2 +2r tf cos (f-0) 
l)—(z„ \-mr 0 cos \p | m cos 0)/ v ; l -j 



NORMAL COMPONENT OF INDUCED VELOCITY FOR A UNIFORMLY LOADED LIFTING ROTOR 


o 



Tlio Wilke anglo * is connected with the resultant 
velocity components at the center of the rotor by 
the relation 

x---tarr , ( mA) (2) 

Jt was desired to compute the nondimensional 
normal eom]>on(*nt of induced velocity in the rotor 
plant 1 at points P(r 0 , \ P) for a wake geometry 
simulating the wake from a rotor operating in the 
upper half of the helicopter speed range. Since 
a wake angle * — tan -1 10 (84.29°) closely approx- 
imates the actual wakt 1 angle for a helicopter 
operating in the higher speed range, the values 
2 o = 0 * 111(1 hi = 1 0 were substituted into equation (1) 
which then became 



2; r 



de 


(3) 


where II represents the integrand of equation (1) 
after the substitutions were made in the quan- 
tities A, B, (\ and D above. 

Numerical approximations to the integral of 
equation (3) were obtained for combinations of 
values of r 0 and ^ by means of a digital computer 
programmed to use Simpson’s rule with 120 equally 
spaced increments in 0. Except for a few points 
close to the wake boundary, this procedure 


yielded results correct to within ±1 in the third 
decimal place as verified by check points previ- 
ously computed by other methods. 

EXPERIMENTAL PROCEDURE 

The electrical systems employed in references 
1 and 2 were broadly similar in that both included 
four basic components : 

(1) The. primary coil (wire model of vortex 

system) 

(2) The secondary coil (search coil) 

(3) The electronic voltmeter 

(4) The power supply 

The methods consisted essentially of measuring 
the voltage induced in the search coil by the 
magnetic field of the primary-coil current and 
converting the result into equivalent velocity. 

In light of information gained from t he reports 
mentioned above 1 , certain fixed considerations 
emerge which affect the accuracy of the method 
and must be taken into account when designing 
an electromagnetic-analog system. These in- 
clude : 

(1) Extraneous magnet ic fields 

(2) Impure wave forms in the primary-coil 

circuit 
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(H) Imluml effects in the primary-coil and 
search-coil leads 

(4) Search-coil dimensions and calibration 

(5) Primary-coil field distortion 

An attorn])! was made in the present, work to 
minimize inaccuracies arising from the above 
sources. The following sections describe each of 
the basic components of the magnetic-analog 
system used in this investigation. 

PRIMARY FIELD ( OIL (WAKE MODEL) 

The difficulties involved in attempting to con- 
struct a solid nonmagnetic cylinder in the shape of 
an elliptic cylinder upon which to wind the pri- 
mary coil made it expedient to build up the wake 
model from a series of “lumped” coils wound on 
separate Plexiglas rings. The rings were mounted 
upon a heavv fiber bast 1 plate by means of indi- 
vidual Plexiglas bases so arranged that the line 
of centers made an angle of 84.29° (or tan 1 10) 
with the rotor plane axis. To minimize the field 
distortion due to lumped coils in the vicinity of 
the rotor plane, the assembly was divided into 
two principal sections. Ihe first section (eone- 
s ponding to the upper portion of the wake) con- 
sisted of 27 rings each bearing 1 turn of No. 17 
gage copper wire. The second section was com- 
prised of IS rings each bearing 9 turns of wire 
and a final ring bearing two layers of 9 turns 
each. The coils were connected in series in 
such a manner that the input and return wires 
for each coil were juxtaposed and could be twisted. 
This arrangement, which for the niultiturn coils 
involved a double w inding, was necessary in order 
to minimize the external magnetic field induced by 
the current in the individual coil leads. The leads 
connecting the wake model to the powci supply 
issued from the final coil at the end of the wake 
model and were also twisted. The wake coils had 
a mean diameter of 12 inches between wire centers 
and were so spaced that the average number of 
turns per unit wake length was the same, in each 
section. It should be noted that the position of 
the rotor plane does not coincide with the plane 
of the end coil but is located approximately half a 
coil turn farther up the wake axis. The relative 
positioning of the coils conformed roughly to the 
actual spacing of the rotor-blade tip vortices in 
the wake of a three-bladed helicopter rotor oper- 
ating at ^.--0.3. The overall length of the assem- 


bly was 12 feet. Under operating conditions the 
“equivalent vortex” strength of the field roil was 
about 4 ampere turns per inch of wake length. 
The entire coil system was mounted on a wooden 
table of height and position such that the wake 
model was centered in its containing room. Figure 
2 is a photograph of the model assembly. 

SEARCH COIL 

The nonlinearity of the primary-coil field and 
the fact that point measurements were desired 
made it necessary that the search-coil dimensions 
be small compart'd with those of the wake model. 
A mean diameter for the search coil amounting to 
about. I percent of that for the field coils was 
adopt el for the work reported herein, since a coil 
of such size could he built with little difficulty and 
would yield induced voltage measurements suffi- 
ciently accurate for engineering purposes. The 
search coil used had a diameter of about 0.35 inch 
to the center of the win* bundle which had a cross 
section in the form of a square approximately 0.09 
inch (>]) a side. The coil consisted of 1.000 turns 
of No. 40 gage copper wire wound on a Plexiglas 
form. The coil form was mounted on a Plexiglas 
suppoit. A solid dielectric coaxial cable was 
list'd to connect the search coil to the amplifier 
in order to minimize the current induced in this 
section of the pickup circuit. The entire search- 
coil assembly together with its coaxial connector 
is shown in figure 2. The base of the search- 
coil support and also the top of the field-coil 
suppoj ting table wen* scribed with straight line's 
spacri at increments of convenient fractions of 
the ro or radius in order to facilitate positioning 
of the search coil. For surveys in the various 
azimuth planes, wooden ramps having the shape 
of 3(T or 60° triangles were used to position 
the seirch-coil assembly. Scribed line's were also 
indue I'd on the faces of these supports. The 
search -coil assembly is shown typically positioned 
relative to the wake model in figure 2. Figure 3 
shows the search coil in detail. 

Tin necessity for obtaining a separate calibra- 
tion e f the search-coil circuit was eliminated in 
the w irk of this report by normalizing the field- 
strong th measurements to those obtained at 
several convenient space locations in the primary- 
coil f old for which the values of the induced 
velocity arc given in reference 4. 
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IwcrKE 2. Magnetic-analog model assembly simulating a wake angle of x tan -1 10 -8-1.29°. 


AMPLIFIER AND OUTPUT METER 

Iii addition to the search coil, the pickup circuit 
included a commercial standing wave indicator 
having a maximum sensitivity of 0.1 microvolt for 
full-scale meter deflection. The assembly con- 
sisted of an indicating meter, a high-gain 400-cycle 
fixed-frequency amplifier with a calibrated gain 
control covering a range of 60 decibels, and a 
narrow 400-cycle band-pass-filter network having 
a sharp cutoff at 400 ± 5 cycles per second. 
The in tegral electronic ally rcgula ted in ten lal 
power supply operated on 115 volts. The input 
impedance of the amplifier was 200,000 ohms and 
consequently it was desirable to test whether 
calibration factors in terms of the search-coil 
current were needed for the indicator readings. 
This was done by placing the search coil at various 


points of high and low field strength and taking 
meter readings with only the normal 200,000- 
ohm impedance in the amplifier input circuit. 
A set of ratios of the equivalent induced velocities 
was computed from these readings. The input 
impedance was then changed to approximately 
5 megohms by means of a noninductive series 
resistor and the procedure was repeated. A com- 
parison of the two siits of computed ratios showed 
no measurable differences. It was concluded 
that meter scale calibration was unnecessary. 
Figure 4 shows the amplifier-indicator unit which 
was located in a hallway removed from the field 
coil. 

POWER SUPPLY 

The power supply used for the wake model 
consisted of a 400-cycle aircraft inverter driven 
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Kk;i hk Details of seairh t-oil. 


by a rectifier, the output voltage of which was 
stabilized by storage batteries. The inverter 
was connected to the primary magnet coil through 
a variable series resistor and through series capac- 
itance. It was found that the frequency stability 
of the system was improved by adjusting the 
capacitance so that the resonant frequency of the 
wake-model coil circuit was slightly above the 
400-cycle operating frequency. As monitoring 
devices the circuit included an ammeter and an 
electrically driven reed frequency meter which 
had been reworked so that the frequencies in- 
dicated by successive reeds differed by only 1 
cps. Rough frequency control was obtained by 
means of inverter taps, and final frequency 
adjustment to the desired 400 cps was made by 
varying the load on the inverter through the 


series resistor. In order to use this frequency 
contrrl system it was necessary to unbalance 
slight lv tin* frequency -load-compensation circuit 
in the inverter. Figure f> shows the power-supply 
assem )ly which was located in a separate room 
from those* of the wake-model coil and amplifier. 

FIELD-SURVEY PROCEDURE 

In general, the wake-model coil circuit was 
allowed to operate for about 30 minutes in order 
to rearh thermal equilibrium before any attempt 
was made to take measurements. After stable 
conditions were reached, the search coil was placed 
at a e.mvenient normalizing point in the magnetic 
field for which the induced velocity ratio was 
known from the digital-computer calculations of 
reference 4 and the meter reading recorded. The 
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Fnii’RK 4. -Fixed-frequency amplifier and indicator unit. 


coil was then moved to the successive survey posi- 
tions and these readings recorded. The search- 
coil circuit was renormalized at frequent time in- 
tervals. 

REDUCTION OF DATA 

The under readings recorded during the pro- 
cedure described in the preceding section were 
converted into equivalent velocity ratios by the 
formula 

! ( * V\ | / 1 Y\ T antilog 0. 1 (AfR) v l 

I \ p 7i> | \ o /.v L an thog 0. 1 (d/AY.vJ 

where Vijr is the nondiinensionul normal compo- 
nent of induced veloi'ity and P refers to tin* spare 
point at which the measurement was made. 

The sign (direction) associated with tin* left 
member of equation (4) was determined from con- 


siderations embracing the flow-field geometry and 
lht‘ trends of the experimental data being reduced. 
The results, as described in the next section, were 
obtained from faired plot s of tin* experimentally 
determined induct'd velocity ratio V f /r plotted 
against R„IR for constant values of ZJR or, where 
necessarv, against Z JR at constant values of 

RalH. 

RESULTS 

Tables 1(a) to 1(g) give the values of 1 \jv 
as experimentally determined over the azimuth 
planes *=0°, 30°, 60°, 90°, 120°, 150°, and 180°. 
Because of the symmetry of the flow, tables 1(b) 
to 1(f) also hold for the azimuth planes ^—330°, 
300°, 270°, 240°, and 210°, respectively. In 
table 1(d) the values of YJv at points for which 


514201 oo 
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Fkurl ">. Power-supply j‘<s« , riihly. 


Og /*>//* g 2.8 and 2 were taken di- 

rectly from the computed results obtained in 
reference 4. The measurements given in table* 
1 extend over a large enough region about tlx* 
rotor plane so that the velocity distributions 
near a rotor with nommiform loading may bo 
calculated by superposition as in reference o. 

Table *2 lists the computed values for V*/t? in the 
rotor plane at azimuth angles \p 0°, 30°, 00°, 
00°, 120°, loO 0 , and 180° extending radially to 
six rotor radii. Although the table contains some 
duplication of values previously listed, it was 
thought convenient in light of possible future 
application to collect the in-plane components 
together. 

Figures 0(a) to 0(d) are plots of constant 
values of Y r ir in tin* various azimuth planes as 


interpolated from t aides 1(a) to 1(g). In par- 
ticular, figure 0(a) supplements the collection of 
similai plots given in reference 3, and figure 0(d) 
extends the ranges covered l>v its corresponding 
plot in reference 4. The* dashed lines in each 
figure represent tin* curve of intersection formed 
by the azimuth plane and the wake vortex cyl- 
inder. Points on these dashed curves are given 
by tin relation 


Xr 


— cot X 




Hr ,, /' (Hr 

/.. eoN*± yi ( /. 


Sin 




(o) 


where only negative values of Z c art* to he. con- 
sidered. 

Figure 7 compares constant -value* plots as ob- 
tained from the computet] values of table 1(d) 
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00 of constant values of Vjr in longitudinal plane 4 , - 0° and 180°. 

itU KK h. Lines of constant values of nondimensional normal component of induced velocity V,7r in each 80° azimuth 
])lane for case of a -wake angle x tan 1 10 8-4.29°. Dashed lines represent cuivh of intersection formed hv azimuth 

plain 1 and wake-vortex cylinder. 
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R 

(,D Li m .s of constant values of T, > in lateral plane 
4 , — 90° and 270°. 

Fir.ru k <». Concluded. 

with those obtained from the experimental values 
in the lateral plant'. The table for these experi- 
mental values has not been included since the 
more accurate computed values were available. 


Figuies 8(a) to 8(e) give plots of the computed 
data of table 2. Experimental analog values for 
the in-plane velocity component are also indicated 
in t best figures for comparison purposes. 

In connection with figures 7 and 8, which show 
compaiisons between computed and experimental 
analog results, it will be noted that no gross dif- 
ferences exist except in regions near the wake 
boundt ry wherein neither the uniform mathe- 
matical model nor t ho magnetic analog with its 
arbitrary finite coil spacing could be expected to 
vield realistic approximations to the true flow 
field. 

CONCLUDING REMARKS 

Inherent in the analog method which has been 
descril ('<1 arc sources of error such as ( 1 ) differences 
in geometry between the model, with its finite 
arbitrary coil spacing, and the wake-vortex system 
for a particular rotor, (2) small variations in 
primary-coil current and frequency, (II) search- 
coil posit ioning errors and associated meter-reading 
errors, (4) inaccuracies in the meter and amplifier 
calibration, and (5) small distortion in tin* portion 
of the model magnetic field of interest arising from 
the laboratory structure. It is to he expected 
that t ic process of fairing the reduced data will 
average out some of the inaccuracies due to the 
above cause's; however, this need not always be 
the ease. Too, the fairing process itself is subject 
to varving degrees of inaccuracy depending upon 
the individual performing the operation. In 
view c f these facts it is difficult to give' any figure' 
for the probable range' of accuracy of the experi- 
mental measurements. However, the' compari- 
sons between the* calculated and analog results 
indicate that the 1 experimental value's are' suf- 
ficiently accurate' for engineering purpose's. 

It is anticipated that the' rompuU'd data pre- 
senter h<T('in will he' useful in synthesizing the' 
(list ributiem of tin* normal component of induced 
ve'loci y ove'r the plane* of any roteer having a, 
specified loading by some method employing tin* 
princ'i )1<* of superposition. Also, it is expected 
that the data should be useful for (estimating the 
int erf* Tcnee'-induee'd velocities of multirotor hcli- 
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(c) Vi/Y in rotor plane for ^- -00° and 150°. 
Fig u re 8 .—Co n eluded. 


rophTs and t ho downwash velocities at wing and 
tail planes. 

Inasmuch as the apparatus and techniques used 
in the present, work are subject to considerable 
refinement, it is thought that the electromagnetic- 
analog method should be useful for mapping 
induced flow fields which are mathematically 
intractable. 

Gkohgia Institute of Technology, 

Atlanta, Ga., February 21 , 1957. 
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NORMAL COMPONENT OF INDUCED VELOCITY FOR A UNIFORMLY LOADED LIFTING ROTOR 
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TAULK 2. NON 1)1 M KN’SION AL VALUES OF NORMAL COMPONENT OF INDUCED VELOCITY U,V IN 
PLANE OF LIFTIN'! I ROTOR FOR WHICH x TAN '1(1 8-1.2!)° 
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